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ABSTRACT 

Using a hybrid method which combines non-radiative hydrodynamical simulations with 
a semi-analytic model of galaxy formation, we determine the ram pressure as a func- 
tion of halocentric distance experienced by galaxies in haloes with virial masses 12.5 ^ 
log (Af 200/1/-^©) < 15.35, for redshifts ^ z ^3. The ram pressure is calculated with 
a self-consistent method which uses the simulation gas particles to obtain the properties of 
the intergalactic medium. The ram pressure profiles obtained can be well described by beta 
profile models, with parameters that depend on redshift and halo virial mass in a simple fash- 
ion. The fitting formulae provided here will prove useful to include ram pressure effects into 
semi-analytic models based on methods which lack gas physics, such as dark matter-only 
simulations or the Press-Schechter formalism. 
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1 INTRODUCTION 

In the local Universe, the properties of galaxies depend on 
the environment in which they reside. The star formation rates 
of galaxies in denser environments such as galaxy clusters are 
lower than for similar field ga laxies |Kauffmann et alj 1 20041 : 
IVerdueo. Ziegler & Ge rken 2008) and the fraction of red galax- 
ies at a given stel l ar mass increases with environm ental density 
jBaldrv et al.l2006l : lMartrnez. Coenda & Muriell 20081) . Disc galax- 
ies in clusters are deficient in HI when compared to similar galaxies 
in the field. The deficiency increases towards the cluster centre, and 
Hi-deficient spirals are o bserved to have truncated gas discs (see 
iBoselli & Gavazzill2006t and references therein). These observa- 
tions suggest that galaxies are transformed from blue star-forming 
systems into red, passive ones via physical processes that remove 
gas, thus suppressing the star formation in the affected galaxies. 

One possible mechanism is ram pressure stripping (RPS) 
of galactic gas caused by the interaction with the hot, high- 
temperature gas of the intragroup or intracluster medium (ICM). 
The galaxy loses gas if the ram pressure (RP) exerted by the ICM, 
P lam = /0icm?J 2 , where v is the velocity of the galaxy relative to 
the ICM, exceeds t he gravitational restoring force of the galaxy 
dGunn & Gottll 19721 : for a recent review of observations and simu- 
lations of RPS in individual galaxies see Roediger 2009). Prior to 
affecting the cold gas discs, RPS could also be responsible for the 
removal of the h ot gas haloes of galax i es after they become satel- 
lites of a cluster ^McCarthy et al J2008I : [Bekkill2009h . 
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Recent results show that RPS is not a process exclusive 
to cluster environments, and it can play a role in smaller sys- 
tems such as galaxy groups (e .g. lRasmussen. Ponman & Mulchaevl 
2006; Rasmussen et al. 2008 1) or for dwarf satellites of giant galax- 
ies (e .g jMc Connachie et al. 2007; Mastropietro, Burkert & Moore 
2009). Thus, a significant number of cluster galaxies could have 
been 'pre-processed' in the small er systems in which they reside d 
prior to infall into the cluster (e.g. Fujita 2004l; ICortese et al.l2006l) . 

Semi-analytic modelling of galaxy formation is a powerful 
approach to explore the effects of environment on galaxy prop- 
erties. Such models do not require heavy computing power, and 
allow the exploration of a large parameter space at a fraction of 
the computational cost of fully self-consistent simulations. Semi- 
analytic models have proven very successful at reproducing several 
observed properties of galaxies such as the local luminosity func- 
tion (see the review by iBaughl 12006). RP S, however, was consid- 
ered i n such models only in a few cases (jOkamoto & Nagashima 
| 2003| : lLanzoni et all 120051 : iBriiggen & De Luciall20o1 " Font et al. 
2008). These models are based on dark matter (DM)-only simula- 
tions, and thus have to resort to analytical approximations to de- 
scribe the ICM. 

An alternative approach is presented by iTecce et alj j2O10l 
hereafter T10), who combine a semi-analytic model with N- 
body /hydrodynamical simulations. In this hybrid approach, the 
kinematical and thermodynamical properties of the ICM are pro- 
vided by the gas particles of the hydrodynamical simulations. This 
results in a self-consistent method which does not introduce ad- 
ditional free parameters into the model, automatically taking into 
account local variations of the density or the velocity field. T10 
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find that, compared to their method, the use of analytic approxima- 
tions results in overestimations of the RP larger than 50 per cent for 
z > 0.5. The disadvantage of the T10 approach is that it cannot be 
combined with DM-only simulation s, such as the widely used Mil- 
lennium simulation (Springel et al. I2005T) . or with m odels which 
use the Extended Press-Schechter (EPS ) formalism dBond et al] 
ll99lUBoweilll99ll ; lLacev & Colelll993l) to determine the growth 
of DM haloes. 

In this paper we determine fitting formulae for the RP exerted 
on galaxies as a function of halocentric distance obtained by using 
the T10 method, in simulated galaxy groups and clusters of differ- 
ent masses. We show that these RP profiles can be well des cribed 
using beta profile models dCavaliere & Fusco-Femianol 19761) . with 
parameters depending on redshift and halo virial mass. The formu- 
lae provided here can then be used in semi-analytic models which 
lack the gas physics required by the T10 approach. 



2 CONSTRUCTION OF THE RAM PRESSURE 
PROFILES 

To calculate the RP exerted by the intergalactic medium on satel- 
lite galaxies in groups and clusters at different redshifts, we use 
the hybrid model sagrp (see T10) which combines non-radiative 
N-body/smoothe d particle hydrodynamics (SPH) simulations of 
galaxy clusters dDolag et alj 1 20051 |2009) with a semi-analytic 
model which handles other baryonic physical processes such as 
gas cooling, star formation, and feedback by supernovae and ac- 
tive galactic nuclei. The SPH simulations are resimulations of the 
regions around five galaxy clusters with mass ~ lO 14 ^- 1 ^ (la- 
beled gl542, g3344, g6212, g676 and g914 in lDolag et all 1 20051 
2009) and three clust ers with mass ~ 10 15 h~ 1 M (g51, gl 
and g8 in iDolag et"aTh . These correspond to a ACDM cosmol- 
ogy with Q m = 0.3, Qa = 0.7, Hubble constant h = 0.7 (in units 
of 100 km s -1 Mpc -1 ), a baryon density fl b = 0.039 and a 
power spectrum normalisation as = 0.9. The DM particle mass is 
1.13 x 10 9 h~ 1 M Q and the gas partic le mass 1.69 x 10 s /i _1 M 
(for more details see IDolag et alj2009h . 

For each DM halo in the sim ulations, id entified by means of a 
friends-of-friends algorithm (FOF. lDavis et alj[l985h , we search for 
all galaxies within the halo virial radius V2oo, defined as the radius 
within which the mean mass density is 200p c ,it- We use these galax- 
ies as tracers of the RP at their current position, calculating for each 
one the RP they experience by using the properties of the surround- 
ing gas particles (see T10 for details). We select all haloes with 
no contamination by boundary particles and with log A/200 ^ 12.5, 
where A/200 is the total mass within r2oo (in units of 1i~ 1 Mq). In 
smaller haloes the RP is very low, affecting only the smallest dwarf 
galaxies, and so these haloes will not be considered in this analysis. 

We identify haloes in the 67 simulation snapshots available in 
the redshift range ^ z ^ 3. The total number of selected haloes 
grows steadily from A^aio = 29 at z = 3 until it reaches a peak 
of 155 haloes at z ~ 0.5, decreasing afterwards to a final value 
of iVhaio = 114 at z = 0. The smallest fof haloes selected are re- 
solved with at least 2800 DM particles, and a similar number of 
gas particles. Haloes with log A/200 — 14 are resolved with about 
9 x 10 4 DM particles, and the most massive haloes are resolved 
with 10 6 DM particles or more. 

Positions and velocities of satellite galaxies within a fof 
halo are determined by the position and velocity of the most 
bound DM particle of the c orresponding subhalo (as done also by 
Briiggen & De Lucia 2008). For those satellites whose subhalo has 
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Figure 1. Top: mean projected galaxy number surface density profile for 
model galaxies with M v < -17 in all haloes with M 2 oo 3 s 10 14 h~ 1 M Q 
at z = 0.3 (solid line). Open circles repres ent the mean o bserved surface 
density profile for CNOC survey clusters ICarlberg et"aT]|l997l) . Bottom: 
mean line-of-sight velocity dispersion profiles for the same haloes, at z = 0. 
Filled and open circles represent, respectively, data f or galaxies with and 
withou t emission lines in the ENACS clusters, from lBiviano & Po ggianti 
<2009h . 



been completely disrupted by the tidal forces exerted by the gravi- 
tational potential of the main halo, we use the most bound particle 
identified at the last time there was a subhalo. This is essentially 
equivalent to assuming that galaxies trace the DM. In the case of 
galaxy clusters, this assumption has been shown to provide a good 
fit to observe d galaxy number profiles out to large cl ustercentric 
radii (see e.g. lBiviano & Girardill2003l:lGao et al1l2004) . 

Using DM particles as tracers for galaxies allows us to gen- 
erate a population of cluster galaxies whose spatial and velocity 
distributions match the observed ones. The top panel in Fig. Q] 
shows the galaxy number surface density profile of simulated 
galaxies with absolute magnitude Mv < —17 in all haloes with 
virial mass ^ 1O 14 /i _1 A/0 at z = 0.3 (solid line), compared 
with the observed profile for clus ter galaxies in the CNOC survey 
dCarlberg. Yee & Ellingsonll 19971 empty circles). The density pro- 
files for the simulations are obtained by projecting along the x, y 
and z axes in turn, considering only galaxies within 2r2oo of the 
cluster centre, binning the galaxies out to 2 r2oo and averaging over 
the three projections (as in lGao et alj|2004l) . The mean galaxy sur- 
face density profile thus obtained agrees very well with the obser- 
vational data. 

The velocity dispersions that we obtain from galaxies con- 
tained within massive haloes (er ~ 400 km s _1 for haloes 
with log Af2oo — 14 and a ~ 900 km s _1 for haloes with 
log A/200 — 15) are consistent with the mass-velocit y dispersion 
scaling relations determined from observations (e.g. lHwang & Lea 
120081: Iwoitak & Lokasll201fj|) . The bottom panel of Fig. □ shows 
the mean line-of-sight velocity dispersion profile of the simulated 
galaxies at z = (solid line), constructed in the same way as the 
density profiles; dashed lines s how the la errors. Circles s how the 
data for ENACS clusters from lBiviano & Poggiantil J2009h . where 
filled and open symbols represent data for galaxies with and with- 
out emission lines, respectively. Although our model galaxies do 
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not present the orbital dicho tomy between early and late type galax- 
ies observed (see e.g. Sodre et alJl989l : lAdarni. Biviano & Mazurd 



1 19981 : IB iviano & Katger tl l2004h . the mean velocity dispersion pro- 
file is a very good match to the average observational trend. There- 
fore, we find that our choice of galaxy orbits is consistent with clus- 
ter observations, and we make the assumption that this method will 
provide a good approximation to the orbits of galaxies in less mas- 
sive haloes as well. 

In T 10 it is shown that the mean RP inside a halo increases 
with redshift and cluster mass. Therefore, we split all the se- 
lected haloes in 10 logarithmic mass intervals in the range 12.5 ^ 
log A/200 < 15.35. The adopted bin width is AlogM2oo = 0.5. 
The bins overlap by half an interval, i.e. 12.5 ^ logM2oo < 13; 
12.75 ^ log A/200 < 13.25 and so on. By doing this we have addi- 
tional points for fitting the mass dependence, in a procedure which 
resembles a (short-period) moving average. For each mass bin we 
construct a combined profile of RP versus halocentric distance by 
gathering the data of all the galaxies within r2oo of the correspond- 
ing haloes. We scale the halocentric galaxy distances by the r2oo of 
their host halo, and we determine the median RP in equally-spaced 
radial intervals. 

To determine whether the shapes of the RP profiles obtained 
depend on virial mass and redshift, we renormalise them to their 
value at xn = rjv/r^oo = 0.15. For the less massive haloes this is 
equivalent to 6 times the gravitational softening used in the SPH 
simulations. As an example, Fig. [2] shows the combined RP pro- 
files determined for four different mass bins at z = 1 (top panel) 
and 2 = (bottom panel). As can be seen, the combined RP pro- 
files show different slopes for different mass intervals and times. 
The RP gradient is stronger in the more massive clusters, where 
the RP at r = T2oo is ~ 1 per cent of the value at xn\ in the 
less massive haloes, the RP at r2oo is ~10 per cent of the cen- 
tral value. We also find, at a given halo mass, a dependence on 
redshift of the normalisation of the RP profiles. For example, the 
average RP value at xn varies from ~1.8 x 10~ 10 h 2 dyn cm -2 
at z = 1 to ~4.7 x 10~ 10 h 2 dyn cm -2 at z = for haloes with 
14.75 «S log A/ 2 oo < 15.35, and from ~2.1 x 10" 12 h 2 dyn cm -2 
at z = 1 to ~4.6 x 10~ 12 h 2 dyn cm -2 at z = for haloes with 
12.5 ^ log A/200 < 13. These trends suggest the need to adopt 
a fitting function with at least two parameters, normalisation and 
slope, both depending on halo virial mass and redshift. 



3 FITS TO THE PROFILES 

To fit the combined RP profiles obtained at each simulation snap- 
shot for the different halo mass ranges selected, three different an- 
alytical functions are considered: a simple power law, 

Pram oc (r/rf)" a , (1) 

a iNavarro. Frenk & White! <1997i hereafter NFW) profile, which 
provides a good description of the density profiles of DM haloes, 

P lm oc (r/rf)- 1 [l + (r/rf )] ~' , (2) 

and a beta profile, which is commonly used to fit the density profiles 
of the ICM in galaxy clusters, 

PramOC [l+(r/rf) 2 ] _3,3/2 . (3) 

In the above relations rf , ry an d r f are the corresponding charac- 
teristic radii. When used to describe ICM profiles, the beta model 
is usually chosen to have a fixed exponent j3 = 2/3. In our case, 
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Figure 2. Combined RP profiles for haloes in four selected virial mass 
ranges at z = 1 (top) and z = (bottom). Each profile is normalised to 
its value at rjv A"200 = 0.15, which for the less massive haloes is equivalent 
to 6 times the gravitational softening used in the N-body/hydrodynamical 
simulations. 



we fit both a model with this fixed value, and another one where f3 
is allowed to vary; these models will be hereafter called fixed beta 
and full beta model, respectively. 

We determine the best-fitting profile of the form (TJl, l[2} 
and ((3) (in this latter case, for both the fixed and the full beta mod- 
els) for the combined RP profiles in the ten mass intervals consid- 
ered and as a function of redshift. All fittings are carried out us- 
ing the L evenberg-Marquar dt technique to solve the least-squares 
problem (Markwardt 2009). For each fit, we compute the chi- 
square goodness-of-fit estimator. The chi-square values are plotted 
in Fig.[3]as a function of the cosmological expansion factor a ex p of 
the corresponding simulation snapshots, for three mass ranges se- 
lected for illustration purposes only. To make the comparison eas- 
ier, the chi-square values in Fig.[3]have all been normalised to the 
average chi-square value obtained for the full beta model, which as 
can be clearly seen from the figure, provides the best-fitting profile 
for all times and mass ranges. Therefore, we propose as a model 
for the RP profiles determined with the T10 method a full beta 
model (O with a core value Po, scale radius r s and exponent f3 
all depending on both virial mass and redshift, 



P mm (M,z) = P (M,z) 



1 + 



,{M,z) 



\f){M,z) 



(4) 



The next step is to determine how the model parameters in ® 
depend on DM halo mass and redshift. To do this we plot for each 
mass interval the best-fitting parameters as a function of a exp . We 
find that for all the mass intervals considered, the evolution of the 
best-fitting parameters can be well fitted with a linear regression. 
As an example, Fig. [4] shows the evolution with time of the best- 
fitting values of logp), r s and f3 obtained for haloes in the range 
13 < log A/200 < 13.5. 
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Figure 3. Chi-square goodness-of-fit estimator as a function of cosmologi- 
cal expansion factor a eX p for the four analytic profile models considered, 
and for three selected mass ranges: 14.75 C logAf 2 oo < 15.35 (top), 
13.75 logAf 20 o < 14.25 (centre) and 12.75 sg logA/ 2 oo < 13.25 (bot- 
tom). In each mass range, x 2 values for all models are scaled to the average 
value obtained for the full beta model. 
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Figure 4. Evolution with time of the best-fitting parameters for the full beta 
model, for the mass range 13 ^ log M200 < 13.5 selected as an example: 
logarithm of the core value Pq (top), scale radius r s /V200 (centre) and ex- 
ponent /3 (bottom). The dashed lines indicate linear fits to the corresponding 
parameters, of the form A + B(a eX p — 0.25). 



The zero point and slope of the linear fits change for the dif- 
ferent mass ranges. Hence, we propose a dependence of the model 
parameters of the following form: 



log 



Po 



10- 12 /i 2 dyn cm" 2 

f"200 



= A P + B P (a cxp ~ 0.25), (5a) 
= A. + B,.(aex P -0.25), (5b) 
= ^ + B^(a cxp -0.25), (5c) 



where in the above relations, the coefficients A and B all depend, 
in principle, on halo virial mass. We have chosen to express the RP 
values in units of 10 -12 h 2 dyn cm -2 , which is of the order of the 
core value of RP in the smallest haloes considered. 

The fitting coefficients for the linear regressions (5} corre- 
sponding to each virial mass range are shown in Fig.[5] In all cases, 
the dependence of the linear fit coefficients A and B on virial mass 
can again be very well fitted by a linear regression, chosen to be 
of the form a + 6(log M200 — 12). Combining the mass and redshift 
dependences, we finally obtain the following expressions for the 
coefficients in ((5): 



A P 


= (-0.8 ±0.1) + (1.2±0.1)(logM 20 o - 12), 




(6a) 


B P 


= (1.2 ±0.2) ± (-0.4±0.1)(logM 20 o - 12), 




(6b) 


A r 


= (0.59 ± 0.03) ± (-0.14 ± 0.02) (log M 20 o - 


12), 


(6c) 


B r 


= (-0.44 ± 0.06) ± (0.12 ± 0.04) (log M 20 o - 


12), 


(6d) 


Ap 


= 0.92 ±0.08, 




(6e) 




= -0.4 ±0.1. 




(6f) 



To illustrate the potential of our model, we have run a modi- 
fied version of sagrp in which we replaced the RP calculation us- 
ing the gas particles from the underlying simulations by the fit- 
ting formulae determined in this work. We will hereafter refer to 
this modified model as sagrp-f. In sagrp-f, as an approximation we 
set the RP equal to zero for all galaxies if z > 3, and for 2 = 3 
onwards we determine, for all satellite galaxies in each halo, the 
distance to the halo central galaxy and then the corresponding RP 
using 10, and ((6). We compare the results of sagrp -f with those 
from sagrp and also with sagrp-a, a model where the RP is calcu- 
lated analytically by using a NFW profile for the ICM density and 
assuming the ICM to be in hydrostatical equilibrium (see T 10 for 
details). 

Fig. [6] shows the fraction of satellite galaxies with stellar 
mass ^ 10 9 h~ 1 M ( r ) that have completely lost their cold gas 
mass as a function of halocentric distance, for the three differ- 
ent models runs: sagrp (solid lines), sagrp-a (dot-dashed lines), 
and sagrp-f (dashed lines). Results are shown grouping the halo 
masses in three separate bins: 14.5 ^ logM 2 oo ^ 15.5, 13.5 ^ 
log Af 2 oo < 14.5 and 12.5 ^ log M 2 oo < 13.5 (top, centre and bot- 
tom rows, respectively), and for three selected redshifts: z = 1 (left 
column), z = 0.5 (centre column) and z = (right column). 

As already shown in T10, calculating the RP by using an an- 
alytical approximation is a good match for the most massive clus- 
ters at 2 ~ 0, but grossly overestimates the RP for log A/ 2 oo < 14 
at 2 > 0. On the other hand, Fig. [6] clearly shows that our best- 
fitting formulae provide an excellent match to the results of the 
self-consistent RP estimation of sagrp, matching the slope and evo- 
lution of the fractions in all mass ranges and at all the redshifts 
considered. For the two most massive ranges considered (top and 
centre rows in Fig.[6]l, the sagrp-f model appears to systematically 
underestimate the RP; however, the diference is within the lcr er- 
rors for the models. 
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Figure 5. Values of the parameters obtained when fitting the time evolution of the best-fitting parameters for the full beta model (4) with a straight line 
A+B(ciz X p— 0.25), plotted as a function of halo virial mass: log Po (top panels), r s /riao (centre) and (bottom). Left and right panels show the corresponding 
values for the coefficients A and B, respectively. In all cases, the mass dependence of the coefficients can be fitted again with a linear regression of the form 
a + 6(log M200 — 12.0), shown with a dashed line. 
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Figure 6. Fraction of galaxies with stellar mass ^ 10 9 /i _1 Mq that have completely lost their cold gas as a function of halocentric distance, for three 
different models: SAGRP, which determines the RP from the information provided by the gas particles in the underlying simulations (solid lines), SAGRP- 
A, which calculates the RP by using an analytical approximation (dot-dashed lines), and SAGRP-F, which uses the best-fitting formulae determined in this 
work (dashed lines). Results are presented for 2 = 1, 0.5 and (left, centre and right columns, respectively) and for haloes in three different mass ranges: 
14.5 log M200 ^ 15.5 (top row), 13.5 ^ log M200 < 14.5 (centre row) and 12.5 ^ log Afjoo < 13.5 (bottom row). Error bars show the la errors for the 
SAGRP model. 
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4 CONCLUSIONS 

Using a hybrid method that combines a semi-analytic model of 
galaxy formation with cosmological hydrodynamical simulations, 
which takes advantage of the extra kinematical and fhermodynam- 
ical information provided by the simulation gas particles, we have 
determined the RP as a function of halocentric distance for DM 
haloes with virial masses 12.5 ^ log M200 < 15.35, for redshifts in 
the range ^ z ^ 3. The RP profiles found can be well described 
by beta models, with parameters that depend on virial mass and 
redshift in a relatively simple fashion. The resulting prescriptions 
can be used to include the effect of RPS in semi-analytic models in 
which the growth of DM haloes is determined either from cosmo- 
logical N-body simulations which do not includ e gas particles, or 
with an approach based in the EPS formalism (e.g.lShefh & Tormenl 
l2002tlCole et alJl2008l ; lNeistein. Maccio & Dekefeoid) . 

In our best-fitting formulae for the RP profiles given by I0, 
(|5j and (|6), we note that whereas the core RP value Po and scale 
radius r 3 depend on both halo virial mass and redshift, the exponent 
/3 depends only on time. We have checked this by running several 
different variants of the sagrp-f model, including in equations < [6et 
and J6fj a mass-dependent term similar to those in the relations 
for the other parameters, varying the coefficients bA « and fes^ (see 
Fig. [5} within the range of errors obtained for them. In all cases, 
we find that the smallest difference between sagrp and sagrp-f is 
obtained when /3 is assumed to be independent of virial mass. 

The results pres e nted i n this work extend those of T10 and 
Briigg en" & De Lucial d2008h . who focused on the most massive 
clusters with log M200 > 14, by also considering the RP profiles 
in less massive galaxy groups. The gradient of RP becomes steeper 
as the virial mass increases; in the massive clusters, the RP at the 
core is ~100 times higher than at r = r2oo, whereas in galaxy 
group-sized haloes the RP at 7-200 is about 10 per cent of the central 
value. On the other hand, the RP in the outskirts of log A/200 =15 
clusters is ~5 x 10 -12 h 2 dyn cm -2 at z = 0, of the same order 
of magnitude as the RP experienced by galaxies in the cores of 
haloes w ith logM2oo = 13.5 a t the s ame epoch. Using SPH sim- 
ulations, Roediger & Briiggen (2006) find that such levels of RP 
can remove about a quarter of the total gas of a spiral galaxy with 
mass ~2 x 10 11 Mq. One concludes that dwarf galaxies with stel- 
lar mass < 10 9 Mq in groups of all masses are very likely to expe- 
rience gas loss by RPS, at least in the local Universe when the RP 
levels are higher in all haloes. 

The fits obtained in this work were determined by using a set 
of resimulations including gas physics of the regions surrounding 
eight massive galaxy clusters, extracted originally from a DM-only 
cosmological simulation. These regions are large enough to contain 
a fair sample of haloes in the mass range considered, free of con- 
tamination from boundary particles. Haloes in the lowest mass bin 
are resolved with at least 2800 DM particles, and the most massive 
haloes are resolved with 10 6 DM particles or more. However, we 
are restricted to only a few massive clusters (there are only three 
haloes in the most massive bin), and this increases the relative error 
in the fitting parameters. 

The RPS depends on the orbits of the satellite galaxies. There 
are several different approaches to dete rmining the orbits of satel- 
lites i n a semi- analyt ic mod el (see e.g. iLan zoni et al. 2005j; ICoral 
120061 : Font et al.|2008). As in lBruggen & De Lucial d2008l) . we track 
the orbital evolution of satellite galaxies within haloes by follow- 
ing the most-bound DM particle of their subhalo. We have shown 
that the method chosen results in positions and velocities of clus- 
ter galaxies which are in good agreement with the mean observed 



trends, although in the case of the velocity dispersion profiles the 
dichotomy between early and late type galaxies is not seen in the 
model galaxies. This may point to missing physics in the models, 
since it is still u nclear which process is resp onsible for the observed 
differences (see lBiviano & P oggianti 2 0091) . 

The particular profile fits obtained could also depend on the 
method chosen for the hydrodynamical calculation in the simu- 
lations. For example, the non-radiative SPH simulations used in 
this work may su ffer from artificia lly su ppressed turbulence (see 
iDolag et al.|[2005l and Agertz et al. 2007). Using simulations car- 
ried out using a grid-based method, for example, could result in 
increased turbulence in the ICM and lead to larger scatter in the 
profiles, as already noted in T10. We intend to explore this issue 
in the future, as larger simulations including different implementa- 
tions for the gas physics become available to us. 

Notwithstanding the above considerations, within the current 
knowledge and taking into account the commonly used hypothe- 
ses, our best-fitting formulae capture remarkably well the effects of 
RP as given by the hybrid approach of T10 on galaxies in differ- 
ent DM halo masses at different redshifts, and provide a significant 
improvement over the analytical approximations used so far to es- 
timate RPS effects. 
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